Asymmetric response of a jammed plastic bead raft 
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Fluctuation-dissipation relations have received significant attention as a potential method for 
defining an effective temperature in nonequilibrium systems. The successful development of an ef- 
fective temperature would be an important step in the application of statistical mechanics principles 
to systems driven far from equilibrium. Many of the systems of interest are sufficiently dense that 
they are close to the jamming transition, a point at which interesting correlations develop. Here we 
study the response function in a driven system of plastic beads as a function of the density in order 
to elucidate the impact of the jamming transition on the use of fluctuation-dissipation relations. 
The focus is on measuring the response function for applied shear stress. We find that even when 
the amplitude of the applied stress leads to a linear response in the strain, the time scale of the 
response is dependent on the direction of the applied stress. 

PACS numbers: 05.70.Ln,61.20.Lc,83.50.Ax 
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The identification of an "effective temperature" for sys- 
tems driven far from equilibrium that is analogous to a 
true thermodynamic temperature would represent a ma- 
jor step toward developing a general theory of nonequi- 
librium behavior. One of the reasons that the concept 
of an "effective temperature" is so attractive is the cen- 
trality of real temperature in thermodynamics and sta- 
tistical physics. At its most basic, temperature defines 
when two systems are in thermodynamic equilibrium and 
determines the direction of heat flow, if any, for systems 
in contact. Equipartition provides a connection between 
temperature and average quantities. Temperature enters 
the relation between thermodynamic quantities, such as 
specific heat and derivatives of free energies. Einstein- 
type relations relate various transport coefficients, such 
as diffusion constants and viscosities, through the tem- 
perature. Finally, in linear response theory, tempera- 
ture enters the relationship between time (or frequency) 
dependent response functions and correlation functions. 
The challenge for studies of nonequilibrium systems is 
the determination of which, if any, of the above uses of 
temperature is meaningful as an effective temperature. 
To be a useful concept, a minimum expectation is that 
multiple definitions of effective temperature agree, and 
that there is an understanding of why some definitions 
agree and others do not. 

The use of effective temperatures in driven systems 
has a relatively diverse history. For example, granular 
flows are often characterized by a "granular" tempera- 
ture that is based on the average kinetic energy of the 
particles. Another approach can be described as a "con- 
figurational" temperature. In these studies, (usually with 
dense granular matter) the "atoms" are grains, and the 
phase space is the set of static mechanical configurations 
in which they can be arranged. By considering all possi- 
ble jammed configurations, entropy as a function of den- 
sity is calculated and used to define an effective tem- 
perature |H 0- Q Another approach has been to focus 



on fluctuations in the system as a means of defining a 
temperature. In this case, the probability distribution of 
an appropriate macroscopic variable, such as the power 
or stress, is studied, and the temperature is related to 
various measures of this distribution. This approach has 
been used in turbulent fluid systems |j. Lj. l| & n d driven 
complex fluids pi |a] • 

In this paper, we focus on definitions of effective tem- 
perature based on the relationship between linear re- 
sponse functions and correlation functions 0, Hfj| . There 
are two basic approaches: static response theory relat- 
ing equal time fluctuations with infinite-time response 
and time-dependent methods. Considering a number of 
different static response relations, it was found that mul- 
tiple definitions of effective temperature converge to the 
same value within the context of an athermal model of 
foam 11]. In contrast, the use of time-dependent rela- 
tions have revealed the impact of multiple time scales 
on the dynamics. In simulations of driven, thermal sys- 
tems, it has been shown that the system exhibits the 
real temperature at short times, but it has a higher ef- 
fective temperature on long time scales |l2j . In this case, 
the static effective temperature and time-dependent ef- 
fective temperature do not (and can not) agree with each 
other. A recent surprising result is that when considering 
different conjugate pairs of variables, the static effective 
temperature determined from one pair can agree with 
the time-dependent effective temperature from a differ- 
ent pair [lflj . The conditions for such agreement are not 
yet understood, but recent simulation work for gra nular 
systems suggests interesting directions of study 
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Experiments probing fluctuation-dissipation measures 
in granular systems have considered both diffusivities (a 
static measurement) 01 an d the time-dependent fluctu- 
ation dissipation relationship The results on dif- 
fusivities is consistent with simulations, at least in the 
sense that effective temperatures from equipartition def- 
initions do not tend to agree with fluctuation-dissipation 



2 



definitions if the inertial terms are not relevant. The 
fluctuation-dissipation work considered a single definition 
of temperature for a range of particle probe properties 
|l6| , but it remains to compare different definitions of ef- 
fective temperature experimentally. Finally, an interest- 
ing system that has been studied is ping-pong balls driven 
by flowing air. In this case , th e fluctuation-dissipation 
relation has been confirmed |17j . 

In this paper, we focus on a definition of effective tem- 
perature using a time-dependent linear response measure 
and consider the behavior as a function of density. The 
time dependent method probes both the time-scale de- 
pendence of the effective temperature and the directional 
dependence of the effective temperature. The central is- 
sue is the behavior of the response function as a function 
of the systems density. As the density increases, one 
approaches the jamming transition. The jamming tran- 
sition is the point at which the system becomes "frozen" 
and can be characterized by a diverging viscosity or the 
development of a static elastic shear modulus for the sys- 
tem Our results suggest that the jamming tran- 
sition is a special point where directional effects become 
important. Surprisingly, it is also the point where time- 
dependent linear response measure is most suggestive of 
a meaningful effective temperature. 

The specific system reported on in this paper consists 
of plastic spheres floating on the surface of water. The 
details of the apparatus and measurement techniques are 
found in Ref. 20, 21] . The spheres are confined to the air- 
water interface and placed between two concentric cylin- 
ders. A 50/50 mixture of spheres with radii of 0.32 cm 
and 0.24 cm is used. The packing fraction p is the to- 
tal cross-sectional area of the spheres at their midplane 
divided by the total area of the trough. The radius of 
the outer cylinder is varied between 6 cm and 11.5 cm to 
vary the packing fraction. As jamming occurs with in- 
creasing packing fraction, it is often useful to consider the 
behavior in terms of the inverse of the packing fraction: 
(\> = l/p. Measurements of the critical packing fraction 
for the jamming of bead rafts have been previously re- 
ported I23. Based on these measurements, we focus on 
the range of packing fractions close to, but below, the 
jamming transition. In this regime, the system exhibits 
the properties of a power-law fluid (viscosity proportional 
to the rate of strain to a power) , and there is no evidence 
of a yield-stress. Approaching the jamming transition 
by varying the packing fraction at a fixed rate of strain, 
the average stress, and by definition the viscosity, diverge 
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We are interested in the behavior of the system with an 
imposed constant rate of strain, or flow. Flow is gener- 
ated by rotating the outer cylinder. In order to measure 
the stress generated by the flow, the inner cylinder is 
suspended on a torsion wire and free to rotate. Angular 
displacements of the inner cylinder are measured. The 
stress on the inner cylinder is determined using a — iripi ■ 



Here k is the torsion constant of the wire and r = 2.7 cm 
is the radius of the inner cylinder. Therefore, angular 
fluctuations of the inner cylinder allow for measurement 
of stress fluctuations. From a time series of the angular 
position of the inner cylinder, we can measure correlation 
functions for either the angular displacement of the inner 
cylinder or the stress generated by the flow. 

In addition to its use as a passive measuring device, an 
external torque can be applied to the inner cylinder, and 
the corresponding angular response measured. This is 
achieved through a coil of wire that is attached to the in- 
ner cylinder and suspended within a second coil, referred 
to as the outer coil. A current in the outer coil generates 
a magnetic field, and by applying a voltage to the inner 
coil to change its magnetic moment, a torque can be ap- 
plied to the inner coil [2(| ■ Because the position of the 
inner cylinder is fluctuating during flow, in order to mea- 
sure the response function a repeated square wave pulse is 
applied to the inner coil, and the response for each cycle 
is averaged. For each measurement, the square pulse is 
applied several hundred times, and the results averaged. 
An example of an average response for a 10 s on/10 s 
off pulse is shown in Fig. la. This represents a period of 
20 s. Figure lb shows the individual response in the case 
of pulses with a period of 50 s. To ensure that the sys- 
tem has equilibrated after each on/off pulse, for the rest 
of the paper all results used pulses with a 40 s period. 
This ensures that turning on the pulse corresponds to ap- 
plying a torque in the direction of the external rotation, 
and because the system has fully equilibrated, turning 
off the pulse is equivalent to applying a torque opposite 
the direction of external rotation. 
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FIG. 1: (a) (color online) Averaged angular response (solid 
black line) to an applied square-wave voltage pulse corre- 
sponding to an applied torque (dashed red line). The angle 
and applied torque have been scaled to appear on the same 
plot, (b) Time series of angular response (solid black line) 
given an applied square-wave voltage pulse that corresponds 
to an applied torque (dashed red line) . The angle and applied 
torque have been scaled to appear on the same plot. 



We will follow the conventions of Ref. [13J and consider 
the correlation function for the angular displacements of 
the inner cylinder, Cg(t) =< 0(i)0(O) > - < > 2 , and 
the integrated response of to a small, constant applied 
torque 5t given by Rg(t) = [< 6{t) > - < 0(0) >]/Sr. 
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Defining rescaled quantities, R$(t) = Rg(t)/Cg(0) and 
Cg{t) = Cg(t)/Cg(0), the fluctuation dissipation relation 
gives Rg(t) = (1/T)(1 - C 6 (t)), where T is the temper- 
ature, or in our case, the effective temperature. It is 
worth briefly considering the case of an overdamped, tor- 
sional oscillator in a thermal bath described by Langevin 
equation: kO + a6 = £(t), where k is the torsion con- 
stant, a is the damping coefficient and £(t) is the ran- 
dom force term for which < £(t)£(0) >= 2aTS(t). For 
such a system, C$(t) = (T/k) exp[— (fc/a)t] and Rg(t) = 
(1/k)(1 — exp[—(n/a)t). From this, we see that the fluc- 
tuation dissipation relation is a statement that the relax- 
ation times for the correlation and response function are 
the same. 
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FIG. 2: The symbols are the final averaged angular displace- 
ment of the inner cylinder as a function of the applied voltage 
pulse given an external rotation rate of 0.01 rad/s and a pack- 
ing fraction of 0.75. The applied torque is linear in the applied 
voltage.) 



To test for linear response, we applied square waves of 
increasingly larger amplitude voltages to the inner coil, 
generating a torque on the inner cylinder. Figure 2 is a 
plot of the final angular displacement of the inner cylin- 
der versus applied voltage for a packing fraction of 0.75 
and a rotation rate of 0.01 rad/s. This illustrates the lin- 
earity of the response. For the range of voltages used, the 
response functions all collapsed onto a single curve when 
scaled by the applied voltage, confirming linear response. 
It should be noted that turning on the voltage pulse cor- 
responded to decreasing the applied rate of strain. Like- 
wise, turning off the pulse corresponded to increasing the 
applied rate of strain. 

Before considering the connection between the re- 
sponse function and the correlation function, it is use- 
ful to consider the time scales associated with each sep- 
arately. Figure 3 is a plot of the relaxation times for 



FIG. 3: (color online) Plot of the response times as a function 
of rotation rate for three different packing fractions. Solid 
symbols correspond to applying a torque in the direction of 
flow, and open symbols correspond to a torque opposite the 
flow. The insert is a plot of the correlation times as a function 
of rotation rate for the same three packing fractions. 



different densities and rotation rates that were obtained 
by fitting the response function to an exponential curve. 
For comparison, the insert shows the correlation times 
for the same densities as a function of rotation rate. It 
should be noted that the correlation function in time was 
well fit by an exponential (in frequency space, it was 
Lorentzian). This behavior is consistent with modelling 
the inner cylinder as an over-damped harmonic oscillator 
in a bath. The main feature of Fig. 3 is the fact that the 
discrepancy between relaxation times for applied torques 
parallel and anti-parallel to the applied rotation rate de- 
velops only for high densities and slow rotation rates. For 
this system, as the rotation rate is decreased, the aver- 
age stress during flow is decreased jl^. The jamming 
transition is approached either by increasing the pack- 
ing fraction or decreasing the stress Therefore, the 
anisotropy in response develops as the jamming transi- 
tion is approached. 

To highlight the features of the response and correla- 
tion functions illustrated in Fig. 3, Fig. 4a consists of a 
series of parametric plots of Rg (t) as a function of Cg (t) 
for packing fractions, 0.72, 0.74, and 0.78, respectively. 
The solid lines are for the response when the applied 
torque is in the direction of flow (decreasing the rate of 
strain). The dashed lines are for the response when the 
applied torque is in the opposite direction as the rotation 
of the outer cylinder (increasing rate of strain). If the sys- 
tem had a single effective temperature, these plots would 
be straight lines with a slope of 1/T and a y-intercept 
of 1/T. Also, if the direction of applied torque did not 
matter, the two curves for each density would superim- 
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FIG. 4: (a) Parametric plots of the angular response function 
versus the angular correlation function. The solid lines are 
for the response when the applied torque is in the direction 
of flow. The dashed lines are for the response when the ap- 
plied torque is in the opposite direction of the flow. In order 
of decreasing amplitude for the response, the curves are for 
packing fractions of 0.72, 0.74, and 0.78. (b) A close-up of the 
two curves for packing fraction 0.78. The dashed and solid 
lines are the same curves as in (a). 



pose. As expected, as the density is increased, the curves 
do not superimpose because the relaxation times are dif- 
ferent for the different directions of applied torques (see 
Fig. 3). 

In terms of defining an effective temperature, it is in- 
teresting to note that as the packing fraction is increased, 
the response as a function of the correlation approaches 
a straight line for the case of a torque applied in the same 
direction as the flow. This is illustrated in Fig. 4b where 
the results for a packing fraction of 0.78 are plotted. The 
dotted line is a straight line guide to the eye. In this 
case, it is interesting to compare the result this gives for 
an effective temperature with other measures of effective 
temperature to determine what, if any, significance there 
is to the linear behavior. This will be the subject of fu- 
ture studies. 

The results reported here for the measured response 
functions during the flow of plastic beads have a number 
of interesting implications for generalizing fluctuation- 
dissipation relations as methods of defining effective tem- 
peratures. First, they highlight the fact that the applied 
rate of strain can break the symmetry of the system and 
result in a directional dependence for the response func- 
tion. The existence of a directional dependence to the re- 
sponse function is consistent with previous measurements 
of the stress in granular systems in response to changes 
in applied stress or rate of strain |23| . These experiments 



focused on the question of flow start up in an opposite 
direction to a recently stopped applied flow. The initial 
applied rate of strain forms stress chains with a direc- 
tional preference. The difference between flow start up 
in the direction and opposite the direction of the initial 
flow can be understood in terms of the development and 
breaking of these stress chains. This behavior occurs in 
the context of a highly nonlinear change in rate of strain. 
In the studies reported here, we observe similar direc- 
tional dependence, but in a regime where the applied 
torque is small enough that the response is linear in the 
applied torque. An important question is the behavior of 
the system as one continues to decrease the driving am- 
plitude. By definition, the truly linear response function 
can not have the directional asymmetry. We are cur- 
rently working on increasing the resolution of the exper- 
iments to determine if the true linear response function 
is experimentally accessible. However, if the linear re- 
sponse regime is not experimentally achievable, this will 
have important implications for any practical definition 
of effective temperature. 

Finally, this asymmetry in the response function will 
not be a generic feature of all definitions of effective tem- 
perature. For example, experimental measurements in 
another granular system probed the fluctuations perpen- 
dicular to the applied rate of strain, where one would 
not expect anisotropic effects |lGj. It will be important 
for future studies of effective temperatures to determine 
when the directional asymmetry will matter and how this 
impacts the generalization of definitions of effective tem- 
perature in such cases. 
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